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5 ■ ABSTRACT 

(N; 

^ , In this paper, the late-time optical photometry of supernova (SN) 1993J in M81 

^ I from Feburary 1995 to Janurary 2003 is presented. The observations were performed 

I in a set of intermediate-band filters that have the advantage of tracing the strength 

^Sj ' variations of some spectral features. SN 1993J was found to fade very slowly at late 

times, declining only by 0.05±0.02 mag 100 d^^ in most of the filters from 2 to nearly 10 

^ ■ yrs after discovery. Our data suggest that the circumstellar interaction provides most 

^N , of the energy to power the late-time optical emission of SN 1993 J. This is manifested by 

1/-^ , several flux peaks seen in the rough spectral energy distributions constructed from the 

Sy ! multicolor light curves. The flux peaks near 6600 A, 5800 A and 4900 A may correspond 

t:;^- ; to the emission lines of Ha, Na I D + He I A5876 and [O III] AA4959, 5007, respectively. 

^^ ' The evolution of these emission lines suggest a power-law SN density model as proposed 

i-G ■ by Chevalier &: Fransson (1994). 

^ • Subject headings: supernova: individual (SN1993J) — techniques: photometric 

-4— » ■ 

> ■ 

• ,-H , 1. Introduction 

j^ ■ SN 1993J was visually discovered on 1993 March 28 UT by Spanish amateur astronomer 

Francisco Garcia-Diaz (Garcia 1993). It occurred in the nearby galaxy M81 (NGC 3031, d = 
3.63 lb 0.31 Mpc, Freedman et al. 1994) and reached a maximum brightness of My=10.8 mag, and 
is the optically brightest SN in the northern Hemisphere since SN 1954A. SN 1993J is classified as 
a "Type lib" SN: its near-maximum spectra are similar to these of SNe II which are characterized 
by strong hydrogen Balmer lines, and its nebular spectra are similar to these of SNe Ib/c that have 
weak hydrogen but strong He I lines (Filippenko, Matheson, & Ho 1993; Swartz et al. 1993). 
Models of SN 1993J (Nomoto et al. 1993; Podsiadlowski et al. 1993; Woosley et al. 1994) suggest 
that the progenitor has lost all but a small amount of its hydrogen layers, thus the shock-heated 
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effective photosphere could quickly sink through the thin H layer into the deeper He layers during 
the initial expansion and the cooling phase. The light curve of SN 1993J was unlike either of the 
plateau or linear version of the typical SN II, showing two maxima similar to SN 1987A but with a 
faster evolution. Except for the initial peak, the light curve bears a strong resemblance to SN Ib/c. 

Due to its brilliant brightness and peculiarity, SN 1993J has been extensively observed along 
the whole electromagnetic spectrum with the telescopes from the ground and in orbit. Early 
observations of photometry and spectroscopy in the optical bands are reported by many authors 
(Schmidt et al. 1993; Wheeler et al. 1993; Filippenko, Matheson, & Barth 1994; Lewis et al. 
1994; Benson et al. 1994; Richmond et al. 1994; Prabhu et al. 1995). The optical photometry 
during the first 360 days after explosion are published by Barbon et al. (1995) and Richmond et 
al. (1996). The optical spectroscopy during SN 1993 J's first 2500 days is presented by Matheson 
et al. (2000a,b). 

Late-time observations of SNe II have established that some SNe II halt luminosity decline, 
remain optically detectable for years and even decade after an outburst. Continued observations 
of these old, evolved SNe can make constraints on the energy mechanism of the late-time optical 
emission and provide clues to the evolution of their progenitor systems. 

In this paper we present the late-time photometric observations of SN 1993J from the years 
1995 to 2003. We describe our observations and data reduction in §2, and the multi-color light 
curves and spectral energy distribution (SED) in §3. Discussions on the late-time energy mechanism 
of SN 1993J are presented in §4, and the results are summarized in §5. 



2. Observations and data reduction 

Photometric observations of SN 1993J at late times have been obtained with the 60/90 cm 
f/3 Schmidt telescope located at the Xinglong station of the National Astronomical Observatory of 
China (NAOC). A Ford Aerospace 2048x2048 CCD camera with a 15 micron pixel size is mounted 
at the Schmidt focus of the telescope. The field of view of the CCD is 58' x 58' with a plate scale 
of 1.67 arcsec per pixel. 

This telescope has a photometric system with 15 intermediate-band (FWHM ~ 200-400 A) 
filters covering a wavelength range from 3000 A to 10000 A (Fan et al. 1996; Yan et al. 1999; Zhou 
et al. 2003). The transmission curves of these filters are shown as the dashed curves in Figure 
1. These filters which were designed to avoid contamination from the strongest and most variable 
night sky emission lines, were used to conduct a survey project among astronomers in Beijing, 
Arizona, Taipei, and Connecticut (BATC) (Fan et al. 1996). We adopt the spectrophotometric 
AB magnitude system in which the fiux calibration is performed by observing four F sub-dwarfs, 
HD 19445, HD 84937, BD-F262606, and BD-M74708 (Oke k Gunn 1983). The transformation 
equations between the broad-band Johnson-Cousins UBVRI and our intermediate-band are the 
following (Zhou et al. 2003): 



U = 6 + 0.6801(a- 6) -0.8982 ±0.143, 

B = d + 0.2201(c-(i)+ 0.1278 ±0.076, 

V = g + 0.3292(/ - h) + 0.0476 ± 0.027, (1) 

R = i± 0.1036 ±0.055, 

I = o±0.7190(n-j9)- 0.2994 ±0.064. 

We have observed SN 1993J roughly once a year using 12 filters {d to p) since 1995. The details 
of the observations are given in Table 1. 



2.1. Photometry 

To obtain proper photometry of an SN that occurs in a complicated background (e.g. spiral 
arms or H II regions), observers are usually required to take template images of the host galaxy 
long after the SN has faded and do image subtraction. This method is appropriate when the SN has 
a relatively fast evolution and the template images can be obtained in a reasonable time. However, 
some SNe, as the case for SN 1993J, are long-lived and evolve slowly during the nebular phase, 
and the template images cannot be obtained without the contamination of the SN light (Li et al. 
2002). Therefore, we have to find an alternative method to fit the galaxy background underneath 
SN 1993J. 

We have used a method that assumes the spiral plane around SN 1993J in M81 is a stable 
system which can describe by diffusive eqution. This means, if the boundary condition for the 
equation is known, we can get the flux distribution of the spiral arm at the position of SN 1993J, 
and the flux of SN 1993J can be obtained by subtracting the flux of galaxy's spiral arms from the 
total flux. 

Before we numerically solve this Laplace equation, we do a Gaussian smooth to reduce the 
noise. We then estimate the flux distribution of M81 in a circle with a 6-pixel radius centered on SN 
1993 J by interpolating the flux from the pixels outside the circle (see the Appendix). This process 
is repeated until convergency. After subtracting the estimated galaxy flux in the circle, only the 
flux of SN 1993J remains. Figure 2 demonstrates this process. 

The final magnitudes of SN 1993J are measured on the subtracted images with standard 
aperture photometry. We use Pipeline II (a program developed to measure the magnitudes of point 
sources in BATC images) that is based on Stetson's DAOPHOT package (Stetson 1987). 



2.2. Calibration 

A total of 20 photometric nights were used to cahbrate 30 local standard stars in the field 
of SN 1993J. For each photometric night, the afore-mentioned four standard stars were observed 
in a range of airmasses, and we derive iteratively the extinction curves and the slight variation of 
the extinction coefficients K + A{UT). The instrumental magnitudes (rriinst) are calibrated by the 
BATC AB magnitude (mBATc) (Zhou et al. 2001) by 

mBATC = minst + [K + A{UT)]x + C. (2) 

where x is the airmass, and C is the zero point of magnitude. 

Table 1 lists the number of calibrations in each filter, while Table 2 lists the final calibrated 
BATC magnitudes and their uncertainties of the 30 local standard stars (see Figure 3 for a finding 
chart of partial local standard stars). For the purpose of clearly showing the SN image, we intercept 
a picture with the field of view of 30' x 30' from original one degree image, therefore only 5 out of 
30 local standard stars are plotted in Figure 3. 

These local standard stars, being used to calibrate, are listed in Table 3. The estimated error of 
each point is a quadrature of the uncertainties in bias and fiat-field correction, aperture photometry, 
and calibrations (Zhou et al. 2003). The main source of the error comes from photon noises and 
uncertainties in the background subtraction. 

We have also converted the BATC magnitudes of SN 1993J into the Johnson-Cousins system 
{V , R, I bands), in an attempt to connect our measurements to the existing early-time broad-band 
photometry. These converted V, R, I magnitudes are listed in the last three columns in Table 
3. We need to point out, however, the transformation equations we have used (also listed at the 
beginning of this section) are derived from observations of normal stars, and may bear relatively 
large uncertainties when apply to the emission-dominated SN 1993J during its nebular phase. 



3. Multicolor light curves 

3.1. Light curves during the period 1995 to 2003 

In the following sections, we adopt 1993 March 27.5 UT as the day of explosion for SN 1993J 
(Lewis et al. 1994). Our measured photometry of SN 1993 J in the m, n, o, and p bands generally 
have low quality due to the low sensitivity of these filters, and will not be discussed hereafter. The 
light curves of SN 1993J in other bands are shown in Figure 4a (d, e, f, and g bands) and Figure 
4b {h, i, j, and k bands) respectively. 

The slow luminosity decline at late times, as listed in Table 3, is evident for SN 1993J. The 
decline rates in various bands are generally found to be 0.05 =b 0.02 100 d~^ during the years from 



1995 to 2003. The light curves in the k and p bands have only a few points, and reliable decline rates 
cannot be determined. The slow luminosity decline of SN 1993J in all the BATC bands suggests 
that there is a persistent energy source powering SN 1993J up to the most recent observations. 

The light curve of i band which centered on Ha, has overall the best quality among all the 
bands, and displays a distinct two-stage evolution: a relatively faster decline rate of 0.09 it 0.01 mag 
100 d~^ from 700 to 2100 days after explosion, and a relatively slower decline rate of 0.02 it 0.02 
mag 100 d~^ from 2100 to 3600 days after explosion. The light curves in all the other bands are 
consistent with a linear decline at all times. 

Figure 5 shows the evolution of SN 1993J in the V , R, I passbands in the past decade. The 
earlier data are collected from Richmond et al. (1996) and Barbon et al. (1995), while those 
after 500 day were converted from our BATC magnitudes. All the light curves exhibit two distinct 
changes, the first at about 50 day after explosion, and the second at about 400 to 500 days after 
explosion. The late-time (>500 day) decline rate of SN 1993 J in VRI is far smaller than that 
between 50 and 360 days after explosion (0.06 ± 0.02mag 100 d'^ vs 1.75±0.03mag 100 d'^). The 
dramatic slow-down after 500 day after may signal a change of physics for the energy source (see 
§4.2 for more discussions). 

We have also constructed i?-band magnitudes from spectrophotometry based on the spectra 
published by Matheson et al. (2000a, b), and the result is shown as solid circles in Figure 5. The 
i?-band data from spectrophotometry is generally consistent with those converted from the BATC 
magnitudes except perhaps for day 881, which is brighter than the neighboring points by nearly 2 
times. We suspect that this might be caused by they's inaccurate flux calibration. 



3.2. SEDs at several epochs 

The late-time SED of SN 1993J can be best studied by spectroscopy such as these done by 
Matheson et al. (2000a, b). Alternatively, An SED could be constructed from the observed fluxes 
in various passbands at the same epoch. Due to the long exposure times to observe SN 1993J in 
the intermediate bands, however, it is impractical to observe all bands in a single night, and our 
definition of the same epoch refers to a reference date it 30 days. This is reasonable, considering 
the very slow evolution of SN 1993J at late times. 

Figure 6 shows the SEDs of SN 1993J at 700, 1000, 1423, and 3245 days after explosion, 
respectively. The prominent feature in the SEDs is the strong emission from the i band, which is 
an order of magnitude brighter than the neighboring h and j bands. The emissions in the e and g 
bands are also apparent. There is an apparent evolution in the shape of the SED: at 700 to 1000 
days after explosion, the SED peaks at the i band, while at 3245 day after explosion, the e band 
has the strongest emission. 

Figure 1 overplots the transmission curves of the BATC system on the spectral sequence of SN 



1993 J from Matheson et al. (2000b), and it can be seen that our intermediate-band filters have the 
advantage of properly covering some of the flat-topped hnes of SN 1993J. Because of M81's small 
heliocentric velocity (—34 km s~^), the e filter covers the emission lines near 4900 A ([OIII] AA4959, 
5007 and H/3), the g filter covers emission lines of He I A5876 and Na I AA5890, 5896, while the i 
filter covers the strong Ha emission. Thus each of the flux peaks seen in the SEDs corresponds to 
one or more strong emission lines in the spectra. 



4. Discussion 

4.1. The late-time emission lines 

We can gain some knowledge about the evolution of some emission lines from our multicolor 
intermediate-band light curves. In theory, measuring the emission-line intensity directly from the 
photometry is difficult when there is no companion knowledge of the evolution of the line profile. 
This problem is somewhat mitigated because our photometric passbands are intermediate in size 
and the transmission curves are essentially flat in each band. Moreover, the late-time spectral 
sequence published by Matheson et al. (2000b) provides some useful information on the evolution 
of the line profiles. The detected emission lines, in order of decreasing strength, are listed below. 



J).. 1.1. Ha emission 

The total line flux is calculated by summing the contributions from each unit wavelength over a 
range determined by the FWHM. We estimate the FWHM of the Hq line profile at different epochs 
from the published spectra Matheson et al. (2000b). the FWHM of the Ha after 2500 day is 
unknown. According to the angular expansion revealed by The Very Long Baseline Interferometry 
(VLBI) observations (Bartel et al. 2002), the expansion velocity of SN 1993 J's radio shell does not 
change significantly after 1600 day, we therefore assume that the Ha line profile does not change 
after 2500 day. The effect of changing line width on the estimation of Ha flux will not exceed 10% 
even if we allow a decrease of 30 A for the FWHM between 2500 and 3600 days. 

Adopting a distance of 3.63 Mpc (Freedman et al. 1994) and E{B — V) = 0.18 mag toward 
SN 1993 J (Richmond et al. 1994), we estimated the reddening corrected Ha luminosity from the 
i-band light curve. Figure 7 shows the results of our calculations. Note that the calculation does 
not consider the contribution by the continuum. As a result, all measurements suffer from relatively 
large uncertainties. We estimate the uncertainty could be as large as 30% if we consider the fluxes 
measured in the m to p bands as the continuum. 

The current Ha emission corresponds to a luminosity of 4.0 x 10^^ erg s^^. In particular, 
the evolution of the Ha emission from SN 1993J bears a strong resemblance to SN 1980K (Fesen, 
Hurford, & Matonick 1995). The implication of a strong persistent Ha emission at late times for 



SN 1993J will be discussed in S4.2. 



J^.1.2. Eel X5876/Na I XX5890, 5896 emissions 

The gf-band light curve traces the emission lines of He I A5876 and/or Na I AA5890, 5896. The 
luminosity of the blended lines fades very slowly from year 1995 to 2002, decreasing by only ~ 1.2 
mag. The total flux measured at the last observation (2002 Mar 3 UT) is 1.8 x 10~^^ erg cm^^ 
s~^, which corresponds to a luminosity of 4.6 x 10^^ erg s~^. Detailed analysis of the spectrum on 
day 976 shows that He I A5876 and Na I AA5890, 5896 contribute approximately equally to the 
broad emission feature around 5800A (Matheson et al. 2000b). Assuming the two components 
have relatively the same intensity at all times, the above luminosity suggests that at the age of 8.9 
yr after the SN outburst, the emission line intensity of He I A5876 or Na I D is about half of the 
Ha. 



4.1.3. [O III] A A 4959, 5007 and Hp emissions 

The e-band light curve traces the evolution of [O HI] AA 4959, 5007 and H/3. Matheson et al. 
(2000b) suggested that the contribution of Up to the total flux peak near 4900A may be as high 
as 30-40% before day 976 (e.g. about 40% at day 670), and decreased significantly in the following 
several years (e.g. 22% at day 1766 and 9% at day 2454). As a result, the late-time e-band light 
curve measures primarily the emission from the [O HI] AA 4959, 5007 doublets. At an earlier stage, 
i.e. 2-3 yrs after explosion, the [O HI] line emission is much weaker than Ha and the He I A5876 + 
Na I D blends. At the later stage, i.e. 8-9 yrs after explosion, however, the [O HI] emission seems 
to gain predominance over the other emission lines so that it is the strongest after the age of 8.9 yr. 
The latest e-band measurement (2002 Mar 3 UT) indicates an [O HI] AA 4959, 5007 luminosity of 
5.7 X 10^^ erg s~^, which is higher than the Ha emission at the same epoch. This is also consistent 
with the spectral evolution reported by Matheson et al. (2000b). 



4.2. The energy sources for the late-time optical emissions 

Few SNe have had optical observations up to an age of 10 yr after explosion. The late-time 
photometry, especially through more than one filter, provides useful information on the underlying 
physics for the lingering light, such as the radioactive decay of long-lived isotopes, interaction with 
the circumstellar medium (CSM), light echoes and delayed optical input by finite recombination 
time. The question is how significant each component it is. 

A slow late-time decline of Ha luminosity is in agreement with the predicted energy input 
from ^^Ti decay (Woosely, Pinto, & Hartman 1989). However, "^^Ti or other long-lived radioactive 



species (^^Co, ^'''Co, and ^■^Na) are unlikely to be important energy sources for SN 1993J's optical 
emission at t ~ 10 yr, as the amount of "^^Ti (half-life ~ 60 yr) needed to produce the current 
optical emissions (i.e. i-band emission) is 3.1 x 10~^ (Dm8i/3.63 Mpc)^M0, far greater than that 
predicted from 13 — 25Mq models (Thielemann, Nomoto, & Hashimoto 1996) or produced in SN 
1987A (Woosely, Pinto, & Hartman 1989). Thus the radioactive mechanism is unable to explain 
the bulk of the optical emission of SN 1993J at the age of 10 yr. This means that some other 
mechanisms must dominate the ionization and excitation of hydrogen. 

If the progenitor of SNe emit materials for an extended period of time prior to exploding as a 
SN, they should be surrounded by a dusty CSM. The light from the SN explosion will scatter off 
this dust and produce a light echo. Since SN 1993J exploded near a spiral arm, it is expected to 
illuminate the interstellar and circumstellar material in the form of light echoes. It is likely that a 
light echo component exists in the light curves of many Type II SNe, the question is how significant 
the echo component is (Roscherr &: Schaefer 2000). Sugerman & Crotts (2002) revealed some 
light echo structures around SN 1993J by analyzing archival Hubble Space Telescope data. They 
derived an observed flux from the echoes as ~ 4.3xl0~^^ ergs cm~^ s~^ A~^ in years 2001, which is 
only about 5.4% of the Ha flux derived from our i-band observations at the same time (7.9 x 10^^^ 
ergs cm~^ s~^ A~^). This suggests that the light echoes cannot account for the slow decay of the 
late-time light curves of SN 1993J. 

The recombination emission might give some contribution to the optical input of supernova at 
late times due to the longer recombination time scales, as Kozma & Fransson (1992) showed for 
SN 1987A. However, the delayed recombination should not be responsible for SN 1993J's late-time 
energy. Comparison of the Ha luminosity from model calculations (Kozma & Fransson 1992) and 
observations of SN 1993J shows that the former is far smaller than the latter, e.g. ~1.5x lO^^erg s^^ 
vs. ~2.8xl0^^ erg s~^ at day 700. Moreover, non-thermal emission from a young pulsar appears 
an equally unlikely late-time energy source. Pulsar photoionization nebulae should produce narrow 
emission lines (~ 1000 km s~^. Chevalier & Fransson 1992), in contrast to the broad emission 
features seen in SN 1993 J's late-time spectra (Matheson et al. 2000b). 

Chevalier & Fransson 1994 (hereafter CF94) studied mass loss before the explosion of a SN and 
suggested that the interaction of the SN ejecta with the circumstellar wind could provide persistent 
energy to power the late-time light curves. In the models CF94 studied, cool, freely expanding SN 
ejecta colliding with CSM from a preexpansion stellar wind. A forward shock propagates into the 
wind, while a reverse shock moves back into the ejecta. The SN ejecta have a fairly steep density 
gradient, leading to a slow reverse shock with emission at far-UV wavelengths (possibly in X-rays, 
with a different gradient). This produces emission from highly ionized species. Absorption by a 
shell formed at the shock boundary can yield low-ionized lines, although these can also originate 
in the ejecta themselves. The fact that all late-time, optically detectable SNe (e.g., SNe 1957D, 
1970G, 1979C, 1980K, 1993J) exhibit strong radio emission is consistent with the idea that the 
main late-time energy comes from the interaction between the expanding SN shell and slow moving 
CSM. Two different models of density profile have been considered for the structure of the wind. 



One is a power law, most applicable to a relatively compact progenitor, while the other uses the 
density structure of a red supergiant (RSG) from stellar evolution models. They make specific 
predictions, including the line intensity ratios and line profiles. 

Table 5 list the observed line flux ratio (relative to Ha) and the predicted values from the 
CF94 model, both for the power-law wind structure and the RSG model. The measured flux 
ratios from our data extend from 1.9 to 9.6 yrs after the explosion of SN 1993J, thus providing a 
comparison with several epochs of the models. The flux ratios are corrected for the reddening of 
E{B -V) = 0.18 mag. 

In general, the power-law density structure model provides better fits to the observed line 
ratios than the RSG wind model. In particular, the decrease in Ha emission from 1.9 to 9.9 yrs is 
well produced (the dashed lines in Figure 7) . Assuming Na I D contributes to half of the intensity of 
the emission at 5800 A, its observed flux matches the model prediction at 2 to 5 yrs after explosion, 
but falls short of the prediction at 10 yr (only about 50%). Part of the reason for this discrepancy 
could be caused by the changing ratio between Na I D and He I A5876 at very late times. The 
power-law density structure model of CF94, for example, predicts the Na I D line would be twice 
strong as He I A5876 at late times. Allowing a relatively large contribution of H/3 to the [O III] 
AA4959, 5007-FH/3 blend before day 1766 (20% - 40%), the observed [O III] AA4959, 5007 flux 
matches the model prediction at 2 to 5 yrs after explosion. At 6.7 yr, the contribution of H/3 to [O 
IIIJ+ H/3 blend is less than 10%. If we assume that all flux at 5000 A are caused by [O III] at 10 
yr, we found that the observed [O III] flux does not deviate significantly from the power-law model 
predictions. The RSG wind model (the last column in Table 5B) generally produces much stronger 
lines of [O III] AA4959, 5000 and much weaker Na I D lines than observed. 



5. Conclusions 

We present intermediate-band photometry of SN 1993J from 692 to 3620 days after discovery, 
greatly extending the coverage of its evolution in the optical bands. The intermediate-band light 
curves show a very slow decline after day 700, fading at 0.05 it 0.02 mag 100 d^^. 

We constructed the SEDs of SN 1993J from the measured flux of the multicolor light curves. 
The SEDs show flux peaks near 4900 A, 5800 A, and 6600 A in order of increasing strength, 
which are associated with the line emissions of [O III] AA4959, 5007, Na I D/He I A5876, and 
Ha, respectively. The [O HI] doublet emission gains dominance over time so that it becomes the 
strongest emission at day 3245. 

Several emission lines seen in the SEDs of SN 1993J and their evolution provide evidences that 
the interaction of the ejecta with the GSM is the primary energy source to power the late-time 
optical emission. We also found that the power-law density profile model of the interaction model 
gives quantitative agreement with the observations. In particular, the line ratio of [O HI] doublet 
and Na I D relative to Ha are well reproduced by the model. 
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In the years to come, it will still be interesting to monitor the evolution of the optical emission 
of SN 1993J. In particular, the strength of [O III] AA4959, 5007 will increase continuously according 
to the circumstellar interaction model. Continued observations of old SNe such as SN 1993J offer 
us an opportunity to study the transition of evolved SNe into supernova remnants. 

We are grateful to Mr. Yang Yanbin for his helpful discussions on reduction of the photometric 
data. Financial support for this work has been provided by the National Science Foundation of 
China (NSFC grant 10303002; 10173003) and National Key Basic Research Science Foundation 
(NKBRSF TG199075402). 



A. Appendix: The numerical solution of the Laplace equation 

The Laplace equation is: 

one needs the appropriate boundary conditions to solve this equation for U{x,y). To solve the 
equation on a digital computer, it is usual to discretize the equation and to work with a finite 
lattice. The finite difference version to the Laplace's equation is obtained by using the well known 
approximations : 

U^j^i-2U^j + Uij+i, (A2) 



dx"^ 
dy"^ 



Ui-ij-2Uij + Ui+ij. (A3) 



Here we use i for a row subscript, j for a column subscript, and U is the pixel value at each 
node of the mesh. Adding these two equations and setting the result equal to zero, yields the 
condition that the value at any mesh point must be equal to the average of its neighbors. This is 
another way of defining what is known technically as a 'harmonic' function. Expressed in symbols, 
using four grid neighbors, this gives the equation: 

Ui,j = |(f/i-i,j + Ui+ij + Uij^i + Uij+i). (A4) 
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Fig. 1. — The matching of BATC intermediate-band filters with the SN 1993J's spectra (the late- 
time spectra are taken from Matheson et al. 2000b). 
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Fig. 2. — Comparison of the i-band image (and mesh), and the residual image (and mesh) after 
the fitting background subtraction. Plotted are (a) The original image of SN1993J; (b) The mesh 
of corresponding zone we see in image (a); (c) The fitting background image of the SN; (d) The 
mesh of corresponding zone we see in image (c); (e) The image of SN 1993 J after subtracting the 
fitting background (d); (/) The mesh of corresponding zone we see in image(e). 
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Fig. 3.— i-band BATC image of the field of SN 1993J (M81), taken on 1995 Dec 25. The field of 
view is 30' x 30'. The five out of 30 local standard stars are marked. 
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Fig. 5. — The VRI light curves of SN 1993J over 10 yrs. An offset has been added to V and I 
data for clarity. Overplots (filled circles) are R-band magnitude measured in the late-time spectra 
(Matheson et al 2000b). 
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Fig. 6.— SEDs of SN 1993J at days 700, 1000, 1423 and 3245 after explosion. 
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Fig. 7. — Tlie i-band light curve showing the development of the Ha luminosity. Plotted are (a) 
the Ha luminosity for SN 1993 J over the period 1995-2003 (filled circles), (b) the predicted Ha 
light curve (solid line) from the power law SN density model of Chevalier & Fransson (1994). 
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Table 1. The details of BATC filters and the observations of SN 1993J. 



Filter name 


Ae//(1) 


FWHM(l) 


Limiting (mag) 


iVofe/ 


Nphot'' 


d 


4540 


332 


20.5 


6 


6 


e 


4925 


374 


20.0 


8 


7 


f 


5267 


344 


20.0 


7 


4 


er 


5790 


289 


20.0 


9 


5 


1 


6074 


308 


20.0 


5 


3 


1 


6656 


491 


20.0 


13 


8 


J 


7057 


238 


20.0 


7 


7 


k 


7546 


192 


19.0 


6 


2 


m 


8023 


255 


19.0 


4 


2 


n 


8484 


167 


19.0 


3 


3 


o 


9182 


247 


19.0 


5 


3 


P 


9739 


275 


18.5 


4 


4 



'^The number of images. 
The number of photometric nights. 



Table 2: Photometry of Comparison Stars 



ID. 


a(2000.0) 


<S(2o0o.o) 


d 




c 




f 




s 




h 




i 


J 


k 


in 


n 









P 


1 


09-51:20.65 


68:58:09.3 


13.9981 


'04) 


13.995( 


05) 


13.7791 


:09) 


13.511( 


09) 


13.2631 


:08) 


13.338(021 


13.267 08; 


1 13.128(08) 


13.169(05; 


1 13.1601 


;08) 


13.0131 


;io 


13.025(21 


2 


09:52:10.73 


69:00:10.1 


15.366( 


10) 


15.2871 


10 


15.0651 


26 


14.808( 


23 


14.5521 


23 


14.561(05) 


14.476(21 


1 14.316(23) 


14.330(12 


1 14.2221 


20 


14.2091 


28 


14.153(26) 


3 


09:52:55.42 


69:01:01.6 


14.8741 


07) 


14.849( 


08) 


14.6461 


17) 


14.560( 


18 


14.3291 


19) 


14.397(05) 


14.328(17 


1 14.152(22) 


14.239(10 


1 14.1651 


19 


14.1161 


22 


13.996(17) 


4 


09:59:31.29 


69:01:57.3 


14.0381 


04) 


14.0221 


05 


13.8111 


09) 


13.5831 


10 


13.3641 


08) 


13.422(02) 


13.424(10 


1 13.195(09) 


13.330(06 


1 13.2441 


09 


13.1661 


11 


13.141(24) 


5 


09:51:37.55 


69:03:52.7 


15.6691 


13) 


15.557( 


12 


15.3641 


31) 


15.038( 


27 


14.7721 


26) 


14.78l(06) 


14.797(24 


1 14.516(28) 


14.513(13 


1 14.4761 


24) 


14.3211 


29 


14.357(28) 


6 


10:00:30.51 


69:07:47.0 


16.313( 


20) 


16.3011 


21) 


16.0301 


56) 


15.584( 


42) 


15.2321 


24) 


15.198(08) 


15.043(34 


1 14.848(41) 


14.890(l9 


1 14.8031 


38) 


14.7141 


40) 


14.543(30) 


7 


09:53:16.90 


69:07:02.6 


15.102( 


08) 


14.965( 


■O8) 


14.8311 


21) 


14.244( 


'16) 


13.9421 


13) 


13.892(03) 


13.770(12 


1 13.532(12) 


13.543(06 


1 13.4471 


10 


13.3541 


13 


13.339(27) 


8 


09:53:41.00 


69:08:06.6 


14.1931 


05) 


14.086( 


05) 


13.8951 


09) 


13.685( 


10) 


13.4311 


09) 


13.470(03) 


13.421(09 


1 13.267(09) 


13.296(05 


1 13.2141 


09) 


13.1771 


11 


13.084(23) 


9 


09:56:44.92 


69:09:00.9 


15.521( 


M 


15.5051 


12 


15.2671 


28) 


15.0101 


28) 


14.7591 


26) 


14.739(06) 


14.706(23 


1 14.514(27) 


14.571(14 


1 14.4841 


28 


14.4061 


30 


14.333(33) 


10 


09:57:15.44 


69:11:29.4 


16.9281 


34) 


16.809( 


34 


16.6051 


47) 


15.928( 


34 


15.6981 


31) 


15.681(12) 


15.564(27 


1 15.331(25) 


15.341(27 


1 15.2121 


28) 


14.9481 


18 


14.97l(l6) 


11 


09:59:20.78 


69:13:19.8 


15.619( 


12) 


15.6521 


13) 


15.4491 


37) 


15.2401 


34) 


15.0931 


37) 


15.183(08) 


15.102(35 


1 14.932(23) 


15.025(22 


1 14.9181 


37 


14.7971 


45) 


14.715(33) 


12 


09:59:03.06 


69:13:53.5 


16.890( 


■33) 


16.8581 


36) 


16.6511 


32) 


16.4601 


28) 


15.9511 


27 


16.074(26) 


16.076(21 


1 15.660(l8) 


15.730(28 


1 15.6411 


24 


15.5191 


22 


15.834(30) 
13.572(29) 


13 


09:53:27.32 


69:13:53.6 


14.9821 


08) 


14.846( 


07) 


14.7121 


19 


14.320( 


'15) 


13.9701 


14) 


14.005(04) 


13.947(12 


1 13.768(15) 


13.826(08 


1 13.7441 


15 


13.6371 


17 


14 


09:55:14.89 


69:19:26.4 


14.223( 


05) 


14.2411 


05) 


14.1021 


11 


13.9741 


12 


13.7791 


13) 


13.836(03) 


13.814(11 


1 13.65l(l4) 


13.742(07 


1 13.7091 


14 


13.6411 


17 


13.432(33) 


15 


09:58:12.20 


69:20:11.3 


14.7441 


07) 


14.750( 


07) 


14.6241 


17 


14.358( 


16 


14.1581 


15) 


14.279(04) 


14.309(17 


1 14.076(18) 


14.154(10 


1 14.1081 


19 


14.0501 


23 


li.nsrhs) 


16 


09:58:35.66 


69:20:47.4 


14.9201 


08) 


14.8991 


08) 


14.7891 


20 


14.6101 


19 


14.3761 


20) 


14.456(05) 


14.413(19 


1 14.217(22) 


14.354(11 


1 14.3321 


23 


14.1571 


26 


14.298(27) 


17 


09:56:16.90 


69:20:43.4 


15.4411 


11) 


15.362( 


11) 


15.1321 


25 


14.954( 


28) 


14.7471 


28) 


14.753(06) 


14.700(24 


1 14.564(29) 


14.542(13 


1 14.5401 


26 


14.4571 


31 


14.284(24) 


18 


09:52:54.07 


69:21:15.9 


14.846( 


07) 


14.7831 


07 


14.6601 


19 


14.5231 


I'' 


14.3111 


18) 


14.400(04) 


14.354(17 


1 14.189(22) 


14.321(10 


1 14.2711 


21 


14.1341 


25 


14.280(28) 


19 


09:56:18.84 


69:23:43.7 


15.5081 


12) 


15.470( 


12 


15.3441 


31 


15.135( 


27 


14.8491 


29) 


14.937(07) 


14.818(25 


1 14.705(38) 


14.769(17 


1 14.7731 


29) 


14.7831 


34 


14.520(32) 


20 


09:52:21.53 


69:29:30.6 


15.478( 


12) 


15.3891 


11 


15.3171 


30) 


14.9891 


24 


14.7941 


28) 


14.880(06) 


14.783(26 


1 14.701(35) 


14.723(16 


1 14.6161 


29 


14.6281 


38 


14.616(39) 


21 


09:52:09.61 


69:29:56.3 


17.3351 


'46) 


17.178( 


43 


17.0831 


35) 


16.537( 


39) 


16.6251 


37) 


16.796(37) 


16.493(31 


1 16.497(33) 


16.593(34 


1 16.5301 


34 


16.2191 


29 


16.798(45) 


22 


09:58:47.33 


69:33:14.1 


15.606( 


12) 


15.6401 


12) 


15.4991 


35 


15.3421 


35) 


15.0391 


34) 


15.160(08) 


15.099(33 


1 14.984(39) 


14.940(l8 


1 14.9421 


42 


14.9571 


50) 


14.767(30) 


23 


09:53:49.38 


69:33:56.8 


15.5911 


!l3) 


15.6171 


'12) 


15.3661 


33 


15.4051 


35) 


15.1221 


39) 


15.291(09) 


15.236(35 


1 15.124(35) 


15.192(26 


1 15.1331 


43) 


15.0231 


25 


15.333(34) 


24 


09:59:56.37 


69:39:19.1 


14.6341 


06) 


14.690( 


'07) 


14.4951 


17 


14.349( 


15) 


14.1611 


17) 


14.265(04) 


14.246(17 


1 14.015(19) 


14.139(10 


1 14.1171 


20) 


14.0721 


22) 


14.028(48) 


25 


09:58:21.48 


69:41:17.9 


15.0421 


08) 


14.9891 


08) 


14.7891 


20 


14.5841 


20) 


14.4071 


21) 


14.477(05) 


14.492(19 


1 14.234(22) 


14.297(11 


1 14.3671 


24 


14.2181 


25 


14.125(50) 


26 


10:00:03.31 


69:42:44.3 


16.4371 


24) 


16.227( 


19) 


15.8941 


50) 


15.214( 


32) 


14.9811 


31) 


14.76l(06) 


14.632(22 


1 14.105(20) 


14.054(09 


1 13.9151 


17 


13.7961 


18 


13.675(36) 


27 


09:56:49.93 


69:43:51.6 


17.1251 


41) 


17.1371 


■11 


16.7361 


35 


16.514( 


34 


16.1641 


37) 


16.093(16) 


15.933(23 


1 15.935(26) 


15.710(34 


1 15.5231 


22 


15.4561 


22 


15.700(35) 


28 


09:50:16.86 


69:42:54.6 


16.0291 


'17) 


16.0141 


18 


15.9001 


31 


15.640( 


20) 


15.4221 


25) 


15.509(09) 


15.486(29 


1 15.283(32) 


15.322(28 


1 15.1991 


20) 


15.2691 


19 


14.967(23) 


29 


09:59:21.70 


69:48:37.8 


15.739( 


13) 


15.6991 


13 


15.5041 


37 


15.253( 


35 


15.0331 


35) 


15.056(07) 


15.017(32 


1 14.843(37) 


14.864(17 


1 14.8731 


38 


14.7521 


39 


14.745(35) 


30 


09:54:34.19 


69:48:36.3 


15.108( 


09) 


15.158( 


09) 


14.9271 


:22) 


14.771( 


'23) 


14.5711 


I24 


14.637(05) 


14.604(23; 


1 14.438(26) 


14.464(13; 


1 14.4371 


I24 


14.438( 


32 


14.470(27) 



Note. — All quantities are magnitudes. Uncertainties in the last two digits are indicated in parentheses. 



Table 3: Photometry of Comparison Stars 



Date" 



6fl2 

708 

757 

982 

1000 

1016 

1050 

1105 

1306 

1358 

1423 

1721 

1780 

2059 

2126 

2436 

2487 

2563 

2755 

2888 

3245 

3504 

3620 



18.66(15) 
18.88(20) 
18.96(09) 

19.35(38) 

19.08(46) 
19.19(10) 



18.24(23) 

18.42(10) 
18.39(08) 



18.49(07) 
19.20(61) 



18.81(30) 
18.64(46) 



19.10(26) 



18.87(16) 
19.22(32) 



19.09(64) 
19.50(83) 



17.69(05) 
17.71(08) 



17.91(04) 
18.05(11) 
17.78(05) 

18.15(15) 
18.50(18) 
18.38(21) 

18.62(15) 

19.11(20) 



19.28(11) 
19.44(41) 
19.44(17) 



18.87(30) 
19.10(44) 



19.11(15) 
19.21(58) 
19.43(54) 



18.84(22) 
18.99(41) 



19.24(31) 



19.38(97) 
19.59(13) 
19.72(24) 



18.28(11) 

18.47(18) 
18.56(19) 



18.83(25) 



18.97(35) 
19.25(41) 



19.02(19) 
19.15(24) 



19.53(14) 19.40(12) 20.05(23) 



18.72(26) 



18.82(29) 



19.43(99) 



19.01(94) 
19.54(63) 



19.47(68) 
20.35(100) 



20.70(99) 



19.16(46) 



19.44(37) 



20.11(82) 



20.36(74) 
19.44(41) 



19.14(37) 
19.19(48) 



19.03(54) 



19.17(43) 
18.85(60) 



20.14(96) 



19.19(68) 



19.28(52) 
18.95(57) 



20.28(95) 20.15(110) 20.83(126) 




K3 

CO 



"Days after explosion, 1993 March UT27.5. 
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Table 4. Late-time luminosity decline rates of SN 1993J in 12 intermediate passbands. 



Band 


Phase''(days) 


Decline rate (mag lOOd ^) 


d 


1358-3245 


0.046 ± 0.014 


e 


708-3504 


0.033 ± 0.007 


f 


692-3245 


0.030 ± 0.006 


g 


692-3245 


0.043 ± 0.006 


1 


692-3245 


0.044 ±0.011 


1 


692-3600 


0.066 ± 0.004 


J 


708-3600 


0.041 ± 0.021 


k 


692-3245 


0.071 ± 0.037 


m 


708-3245 


0.048 ± 0.031 


n 


1000-3245 


0.050 ± 0.045 


o 


708-3245 


0.054 ± 0.046 


P 


708-3245 


0.067 ±0.056 



*Days after explosion, 1993 March UT27.5. 



Table 5A. Observed line flux from SN 1993J. 









Age (year) 






Emission Line 


1.9 


2.7 


4.8 5.7 


8.9 


9.6 


TT a 

Na I D±He 15876 
[OI//]4959-5007±i7/3 


2.84 
0.60 
0.58 


2.10 
0.71 
0.69 


1.03 0.63 

0.86 1.12 

••• 1.50 


0.45 
1.06 
1.69 


0.41 

1.72 



Table 5B. Observed line flux from SN 1993J. 







Age (year) 




Emission Line 


2 


5 10 


lO(RSG)'^ 


NalD 
[O///]4959-5007 


0.96 
0.38 
0.27 


0.19 0.09 
0.61 1.00 
1.00 2.20 


0.49 
0.17 
3.40 



^The luminosity of Ha in 10^^ erg s ^. 
'^RSG is the red supergiant model at 10 yr. 



